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Abstract:  The  ADaptive  Hydraulics  model,  ADH,  was  used  to  investigate 
the  circulation  tendencies  in  and  around  Bush  Canal  by  utilizing  its  two- 
dimensional  shallow-water  module.  This  study  characterizes  existing  wa¬ 
ter  levels  and  currents  in  the  vicinity  of  Bush  Canal  and  predicts  any  po¬ 
tential  impacts  that  may  result  from  the  construction  of  a  structure  in 
Bush  Canal  near  its  junction  with  Bayou  Terrebonne.  Comparing  model¬ 
generated  currents  and  water-surface  elevations  between  pre-  and  post¬ 
construction  conditions  provides  insight  into  whether  a  particular  alterna¬ 
tive  will  adversely  impact  velocity  conditions.  ADH  was  used  to  develop 
time  series  of  current  fields  for  the  base  and  plan  conditions. 
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1  Introduction 

1.1  Background 

The  U.S.  Army  Engineer  District,  New  Orleans  (MVN)  is  in  the  process  of 
determining  the  correct  sizing  for  sail-through  structures  to  be  constructed 
along  the  new  proposed  Morganza  to  the  Gulf  levee  system  in  Louisiana. 

In  the  aftermath  of  Hurricane  Katrina,  numerous  new  levees  are  being 
constructed  to  reduce  the  risk  of  storm  surge  flooding.  To  provide  the  best 
protection,  waterways  north  of  the  new  levee  system  will  be  shut  off  from 
the  Gulf  of  Mexico  during  storm  events.  In  order  to  accomplish  this  pro¬ 
tection  with  minimum  impact  on  the  natural  environment  and  waterborne 
transportation,  numerous  sail -through  structures  —  each  capable  of  being 
closed  as  necessary  —  are  being  constructed  to  allow  for  continued  water¬ 
borne  transportation  between  the  Gulf  of  Mexico  and  areas  north  of  the 
new  levee  system.  These  structures  will  also  allow  for  the  continued 
movement  of  water  to  and  from  the  Gulf  of  Mexico,  thereby  minimizing 
the  impact  of  the  new  levees  on  natural  habitat  along  the  Louisiana  coast. 

At  the  request  of  New  Orleans  District,  the  U.S.  Army  Engineer  Research 
and  Development  Center  (ERDC)  performed  a  number  of  engineering 
studies  in  support  of  efforts  to  determine  the  proper  size  of  one  of  these 
proposed  structures,  to  be  located  on  the  eastern  end  of  Bush  Canal  near 
Bayou  Terrebonne.  Four  additional  structures  (Humble  Canal,  Bayou  Ter¬ 
rebonne,  Lapeyrouse  Canal,  and  Bayou  Petit  Caillou)  are  also  to  be  con¬ 
structed  in  the  same  vicinity.  The  Bush  Canal  structure  will  consist  of  a 
sector  gate  structure,  with  additional  sluice  gates  being  added  if  needed. 
The  proposed  levee  alignment  with  all  five  proposed  structure  locations  is 
shown  in  Figure  l. 

1.2  Objective 

The  primary  objective  of  this  study  was  to  determine  the  smallest  structure 
that  will  provide  reasonable  velocity  fields.  The  modeling  approach  ap¬ 
plied  to  size  the  Bush  Canal  structure  will  subsequently  be  used  to  size  the 
other  sail-through  structures. 
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1.3  Approach 

Data  sources  and  acquisition,  and  model  development  methods  using 
ADaptive  Hydraulics  numerical  model  code  (ADH),  are  described  in  detail 
in  Chapter  2.  Several  design  alternatives  were  simulated  using  the  vali¬ 
dated  numerical  model. 


Bayou  Terrebonne  $ 


Canal  Structure 


Bush  Canal  Structure 


Lapeyrouse  Canal  Structure 


Caillou  Structure 


US  SURVEY  MILES 

- @ - 18- 


Figure  1.  The  proposed  levee  alignment  (red  line)  with  the  proposed  structure  locations 
(yellow  circles).  Source:  Google  Maps. 
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2  Model  Development 

2.1  Description  of  site 

The  project  is  located  south  of  Houma,  LA,  on  the  southern  coast  of  Lou¬ 
isiana.  The  primary  waterways  of  interest  are  Bayou  Petit  Caillou,  Bayou 
Terrebonne,  and  Bush  Canal.  Bayou  Petit  Caillou  and  Bayou  Terrebonne 
lie  almost  parallel  to  each  other,  with  approximate  north-south  orienta¬ 
tions.  The  primary  means  of  waterborne  transportation  connecting  Bayou 
Petit  Caillou  and  Bayou  Terrebonne  is  Bush  Canal.  Bayou  Petit  Caillou  ex¬ 
tends  south  to  Cocodrie,  where  it  joins  the  Houma  Navigational  Canal.  It 
also  has  connections  to  the  west  to  Lake  Boudreaux  and  Lake  Quitman. 
Bayou  Terrebonne  extends  south  to  Lake  Barre  and  has  connections  to  the 
east  to  Madison  Bay.  Bayou  Petit  Caillou  and  Bayou  Terrebonne  join  north 
of  the  study  area  and  extend  northward  to  join  the  Gulf  Intracoastal  Wa¬ 
terway  (GIWW)  in  Houma,  LA.  The  study  area  is  shown  in  Figure  2. 
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Figure  2.  Study  area.  Source:  Google  Maps. 
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2.2  Field  data 

2.2.1  Water  surface  elevation  data 

The  first  field  data  collection  goal  for  this  study  was  to  obtain  enough  wa¬ 
ter  surface  elevation  gage  data  to  adequately  validate  the  model  water  sur¬ 
face  elevation  results.  The  validation  period  for  the  model  was  l  May  2006 
to  22  July  2006.  MVN  maintained  several  water  surface  elevation  gages 
located  in  and  around  the  study  area,  and  varying  amounts  of  data  were 
available  for  different  gages  during  this  validation  time  period  (i.e.,  some 
gages  recorded  data  for  the  entire  time  period  and  some  recorded  data  for 
a  shorter  duration  during  this  period).  These  water  surface  elevation  gages 
were  essential  to  the  validation  of  the  model.  The  gage  locations  are  shown 
in  Figure  3. 


Figure  3.  Location  of  the  water  surface  elevation  gages  Source:  Google  Maps. 
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From  inspection  of  the  water  surface  elevation  data  it  was  determined  that 
the  wind  had  a  significant  effect  on  the  water  surface  elevations.  To  re¬ 
move  this  effect,  the  raw  field  data  were  filtered  to  eliminate  any  signals 
that  had  a  period  longer  than  30  hours.  Thus,  the  longer-period  signals 
(primarily  the  wind  signal)  were  removed,  leaving  only  the  shorter-period 
signals  (primarily  the  tidal  signal).  Comparisons  of  the  raw  data  versus  fil¬ 
tered  data  for  each  of  the  gages  are  shown  in  Figures  4  -  10. 


Figure  4.  Graph  of  the  raw  and  filtered  field  data  from  the  Cocodrie  gage. 
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Figure  5.  Graph  of  the  raw  and  filtered  data  from  the  Bush  Canal  gage. 


Figure  6.  Graph  of  the  raw  and  filtered  data  from  the  Lake  Boudreaux  gage. 
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Bayou  Petit  Caillou  Field  Data 
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Figure  7.  Graph  of  the  raw  and  filtered  data  from  the  Bayou  Petit  Caillou  gage. 


Figure  8.  Graph  of  the  raw  and  filtered  data  from  the  Bayou  Terrebonne  gage. 
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Figure  9.  Graph  of  the  raw  and  filtered  data  from  the  Madison  Bay  gage. 


Figure  10.  Graph  of  the  raw  and  filtered  data  from  the  Humble  Canal  gage. 
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2.2.2  Velocity  data 

The  second  data  collection  goal  of  this  study  was  met  by  installing  five  ve¬ 
locity  meters  to  operate  over  the  same  time  period  that  the  surface  eleva¬ 
tion  data  were  collected.  The  purpose  of  the  velocity  meters,  located  as 
shown  in  Figure  n,  was  to  get  a  better  understanding  of  the  flow  split 
along  these  waterways  and  provide  current  velocity  data  for  model  valida¬ 
tion.  The  two  meters  located  to  the  far  north  on  the  two  bayous  (Upper 
Bayou  Terrebonne  and  Upper  Bayou  Petit  Caillou)  were  needed  to  verily 
that  the  tide  was  in  fact  reaching  these  locations  and  that  no  unknown  in¬ 
flows  were  affecting  the  study  area.  These  velocity  meter  data  are  plotted 
in  the  model  validation  section  (section  2.5,  Figures  32  to  34). 


Figure  11.  Locations  of  the  velocity  meters.  Source:  Google  Maps. 
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In  addition  to  the  velocity  measurements,  Acoustic  Doppler  Current  Pro¬ 
file  (ADCP)  transect  measurements  were  taken  on  13  July  2006.  The 
ADCP  measurements  were  taken  at  least  once  an  hour  for  eight  consecu¬ 
tive  hours.  The  locations  for  these  transects  were  chosen  to  allow  for  a 
good  understanding  of  the  flow  splits  occurring  in  and  around  Bush  Canal. 
These  ADCP  measurements  were  used  as  the  primary  validation  tool  for 
comparison  with  the  model  velocities.  The  ADCP  measurements  were  also 
used  to  correlate  the  velocity  meter  measurements  from  Bush  Canal  and 
Bayou  Terrebonne  to  an  average  cross-sectional  velocity.  The  locations  of 
the  ADCP  transects  are  shown  in  Figure  12,  and  the  model  data  are  shown 
in  the  model  validation  section  (section  2.5,  Figures  28  to  31). 


Figure  12.  Locations  of  the  ADCP  transects.  Source:  Google  Maps. 
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The  adjustments  made  to  the  velocities  at  Bush  Canal  and  Bayou  Terre¬ 
bonne  are  shown  in  Figures  13  and  14.  The  velocity  meter  measurements 
at  Bush  Canal  were  increased  by  a  factor  of  2.2  while  the  velocity  meter 
measurements  at  Bayou  Terrebonne  were  increased  by  a  factor  of  3.  These 
factors  were  determined  by  correlating  the  velocity  meter  measurements 
with  the  cross-sectionally  averaged  ADCP  velocities.  By  scaling  the  velocity 
meter  measurements  in  this  way,  a  larger  data  set  is  available  for  model- 
data  comparison.  The  ADCP  measurements  were  relied  upon  in  this  way 
because  they  are  more  representative  of  the  overall  flow  in  the  canal. 


Figure  13.  Plot  of  adjusted  and  unadjusted  velocity  meter  measurements  for  Bush  Canal. 
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Bayou  Terrebonne  ADCP  Velocity  Comparison 


Time  (hrs.  from  Jan  1  0:00,  2006) 


Figure  14.  Plot  of  adjusted  and  unadjusted  velocity  meter  measurements  for  Bush  Canal. 

2.3  Model  mesh  development 

The  ADaptive  Hydraulics  numerical  model  code  (ADH)  was  chosen  to 
model  the  Bush  Canal  system  (see  Appendix).  An  existing  RMA2  mesh  of 
the  central  southern  Louisiana  coast  was  provided  by  New  Orleans  Dis¬ 
trict.  This  initial  mesh  was  created  by  Mr.  Joseph  V.  Letter,  Jr.,  for  an 
Atchafalaya  Bay  study  using  RMA2  (Donnell,  Letter,  and  Teeter  1991).  It 
was  later  modified  by  Mr.  David  Elmore  and  again  by  Amena  Henville 
(both  of  New  Orleans  District)  for  the  Morganza  to  the  Gulf  of  Mexico  pro¬ 
ject,  also  using  RMA2.  This  initial  mesh,  shown  in  Figure  15,  extends  from 
the  Atchafalaya  Bay  to  the  west  to  Port  Fourchon  to  the  east.  It  contains  a 
large  area  to  the  west  of  the  study  area  that  was  not  necessary  for  conduct¬ 
ing  the  current  study. 
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Figure  15.  Initial  RMA2  mesh  developed  for  previous  study,  with  current  study  area  indicated. 

To  create  the  Bush  Canal  mesh  for  the  current  project,  the  western  area 
outside  the  greater  study  location  was  removed.  The  mesh  was  also  con¬ 
verted  to  an  ADH-compatible  format,  which  involved  converting  from 
quadratic  elements  used  in  RMA2  to  linear  elements  for  ADH.  Because 
ADH  accepts  only  triangles,  all  quadrilaterals  in  the  mesh  were  converted 
to  triangles.  The  resolution  in  the  study  area  was  significantly  increased. 
These  mesh  modifications  were  performed  in  the  Surface-water  Modeling 
System  (SMS)  (Brigham  Young  University  2002),  a  graphical  user  inter¬ 
face  developed  by  ERDC  for  use  in  setting  up  and  running  numerical  mod¬ 
els.  The  final  model  domain  used  in  this  study  is  shown  in  Figure  16,  with 
Figure  17  showing  the  increased  resolution  in  the  study  area. 
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Figure  16.  Model  domain  for  current  study. 


Figure  17.  Illustration  of  increased  resolution  in  the  study  area. 
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2.4  Boundary  conditions  development 

Because  there  were  no  river  inflows,  the  development  of  boundary  condi¬ 
tions  consisted  only  of  the  tidal  boundary.  The  predicted  tide  for  Port 
Fourchon  was  applied  along  the  southern  boundary  of  the  model  shown  in 
Figure  18.  Wind  and  rain  were  not  considered  in  this  study,  and  therefore 
were  not  included  in  the  boundary  conditions  file.  The  black  circle  in  the 
lower-right  side  of  Figure  18  shows  the  location  of  the  Port  Fourchon  gage. 


Mesh  Elevation,  ft,  NAVD88(2004.65) 
I  0.0 


Port 

Fourchon 


0.0 


10.0 


20 


.0  3VIILE^0.0  40.0 
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Figure  18.  Mesh  showing  the  tidal  boundary. 


To  ensure  the  validity  of  applying  the  Port  Fourchon  tide  along  the  bound¬ 
ary  farther  to  the  south,  comparisons  were  made  between  the  model  data 
at  the  Port  Fourchon  gage  location  and  the  Port  Fourchon  predicted  tide. 
This  was  done  to  ensure  no  amplification  or  phase  shift  was  occurring  due 
to  the  movement  of  this  tide  south  to  the  model  boundary  location.  A 
graphical  comparison  (Figure  19)  of  the  model  results  at  Port  Fourchon 
and  the  predicted  tide  at  Port  Fourchon  showed  no  significant  difference. 
Therefore,  no  significant  amplification  or  phase  shift  occurs  by  applying 
the  Port  Fourchon  predicted  tide  along  the  southern  boundary  of  the 
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model.  This  is  the  expected  result  since  the  boundary  and  the  gage  are  in 
such  close  proximity  to  each  other.  The  predicted  Port  Fourchon  tide  was 
used  for  the  tidal  boundary  in  the  base  and  plan  condition  model  runs. 


Figure  19.  Comparison  of  model  results  and  predicted  tide  at  Port  Fourchon. 


2.5  Model  validation 

2.5.1  Water  surface  elevation  validation 

The  hydrodynamic  model  was  validated  by  comparison  of  the  model- 
computed  water  surface  elevations  with  filtered  observed  water  surface 
elevations.  New  Orleans  District  had  placed  water  surface  elevation  gages 
in  and  around  the  study  area,  as  previously  shown  in  Figure  3.  These  gages 
were  of  vital  importance  in  the  validation  process.  Comparisons  of  the  fil¬ 
tered  observed  water  surface  elevations  with  the  model  water  surface  ele¬ 
vations  are  shown  in  Figures  20  -  27. 
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From  observing  the  filtered  Cocodrie  data  in  Figure  20,  it  was  noticed  that 
the  range  of  the  tide  changed  dramatically  around  hour  4380.  Upon  fur¬ 
ther  investigation,  it  was  discovered  that  the  Cocodrie  gage  was 
downloaded  at  hour  4382.  After  comparing  model  results  and  the  tides  at 
other  locations,  it  was  determined  that  the  data  prior  to  hour  4382  was 
amplified  by  some  amount.  The  amplified  data  were  then  proportioned 
until  approximately  the  same  spring  tide  range  was  present  for  the  data 
before  and  after  the  gage  reading.  This  proportioning  resulted  in  a  reduc¬ 
tion  in  the  tide  range  (before  hour  4380)  by  a  factor  of  2.  By  halving  the 
range  about  the  mean,  the  tide  ranges  prior  to  hour  4382  approximately 
matched  the  tide  ranges  present  after  hour  4382  for  the  spring  tide  events. 
A  plot  of  the  adjusted  Cocodrie  filtered  data  as  compared  with  the  model 
results  is  shown  in  Figure  21. 
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Cocodrie  Field  Data 


- Cocodrie  Adjusted  Filtered  Data 

- Cocodrie  Model  Results 
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Time  (hrs.  from  Jan  1  0:00.  2006) 


Figure  21.  Adjusted  filtered  data  and  the  model  data  for  Cocodrie. 


Bush  Canal  Field  Data 


Time  (hrs.  from  Jan  1  0:00.  2006) 


Figure  22.  Filtered  data  and  the  model  data  for  Bush  Canal. 
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Figure  23.  Filtered  data  and  the  model  data  for  Lake  Boudreaux. 


Figure  24.  Filtered  data  and  the  model  data  for  Bayou  Petit  Caillou. 
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Figure  25.  Filtered  data  and  the  model  data  for  Bayou  Terrebonne. 


Figure  26.  Filtered  data  and  the  model  data  for  Madison  Bay. 
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With  the  exception  of  the  Lake  Boudreaux  gage,  the  remaining  gage  com¬ 
parisons  are  good.  The  Lake  Boudreaux  gage  is  located  in  Lake  Boudreaux 
near  a  connection  to  Bayou  Petit  Caillou  (see  Figure  3).  Upon  completion 
of  the  hydrodynamic  modeling  it  was  discovered  that  a  weir,  not  included 
in  the  ADH  model,  was  located  along  the  Boudreaux  Canal  connecting 
Bayou  Petit  Caillou  to  Lake  Boudreaux.  This  missing  weir  would  greatly 
affect  the  model  water  surface  elevations  at  the  Lake  Boudreaux  gage.  It  is 
believed,  however,  that  this  missing  weir  will  have  a  minimum  affect  on 
the  remainder  of  the  system  because  the  Bayou  Petit  Caillou  observed  wa¬ 
ter  surface  elevation  and  velocity  data  (Figures  24,  30,  31,  and  32) 
matched  the  model  results  well. 

2.5.2  Velocity  validation 

Validation  of  the  model  velocity  results  consisted  of  comparisons  with  the 
measured  ADCP  data  taken  on  13  July  2006.  The  cross-sectional  averaged 
velocity  in  the  channel  measured  by  the  ADCP  was  compared  to  the  cross- 
sectional  averaged  velocity  in  the  channel  from  the  2D  ADH  model  results. 
The  locations  of  the  ADCP  transects  are  shown  in  Figure  12.  The  ADCP 
comparison  results  are  shown  in  Figures  28  -  31.  The  velocity  meter  loca¬ 
tions  are  shown  in  Figure  11,  and  their  comparison  plots  are  shown  in  Fig¬ 
ures  32  -  34. 
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Bayou  Terrebonne  ADCP  Velocity  Comparison 


Time  (hrs.  from  Jan  1  0:00,  2006) 


Figure  28.  Model  velocity  results  and  the  ADCP  measurements  for  Bayou  Terrebonne. 


Bush  Canal  ADCP  Velocity  Comparison 


Time  (hrs.  from  Jan  1  0:00,  2006) 


Figure  29.  Model  velocity  results  and  the  ADCP  measurements  for  Bush  Canal. 
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Bayou  Petit  Caillou  ADCP  Velocity  Comparison  South  of  Bush  Canal 


Time  (hrs.  from  Jan  1  0:00,  2006) 


Figure  30.  Model  velocity  results  and  the  ADCP  measurements 
for  Bayou  Petit  Caillou  south  of  Bush  Canal. 


Bayou  Petit  Caillou  ADCP  Velocity  Comparison  North  of  Boudreaux  Canal 


Time  (hrs.  from  Jan  1  0:00,  2006) 


Figure  31.  Model  velocity  results  and  the  ADCP  measurements 
for  Bayou  Petit  Caillou  North  of  Boudreaux  Canal. 
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Lower  Bayou  Petit  Caillou  Velocity  Meter  Comparison 


•Velocity  Meter  Measurements  (-  upstream) 
-Model  Velocity  Results  (-  upstream) 
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Time  (hrs.  from  Jan  1  0:00,  2006) 

Figure  32.  Model  velocity  results  and  the  velocity  meter  measurements 
for  lower  Bayou  Petit  Caillou. 


Bush  Canal  Velocity  Meter  Comparison 


Figure  33.  Model  velocity  results  and  the  adjusted  velocity  meter  measurements 

for  Bush  Canal. 
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Bayou  Terrebonne  Velocity  Meter  Comparison 


-Adjusted  Velocity  Meter  Measurements  (+  upstream) 
-Model  Velocity  Results  (+  upstream) 


4200  4300  4400  4500 

Time  (hrs.  from  Jan  1  0:00,  2006) 


Figure  34.  Model  velocity  results  and  the  adjusted  velocity  meter  measurements 

for  Bayou  Terrebonne. 


2.6  Design  alternatives 

The  aim  of  the  design  alternatives  was  to  determine  an  appropriate  struc¬ 
ture  that  would  produce  reasonable  velocity  fields.  Testing  included  the 
use  of  multiple  structure  sizes  in  conjunction  with  the  possible  installation 
of  sluice  gates  of  varying  sizes.  It  was  determined  during  initial  testing  that 
the  design  of  the  structure  and  sluice  gate  entrances  and  exits  had  a  critical 
impact  the  velocities  reported  by  the  model.  By  making  the  sail-through 
structure  more  rounded  at  the  entrance  and  exit,  the  structure  became 
more  efficient  and  passed  more  water.  A  negative  consequence  of  this  de¬ 
sign,  however  was  that  it  also  increased  velocities  within  the  structure.  It 
was  decided  that  the  sluice  gates  should  be  made  more  efficient  by  using 
both  rounded  entrances  and  exits,  but  the  sail-through  structure  should  be 
made  less  efficient  by  incorporating  less-rounded  entrances.  This  design 
alternative  allows  a  significant  portion  of  the  water  to  circumvent  the  less- 
efficient  structure  and  flow  instead  through  the  more-efficient  sluice  gates. 
This  results  in  increased  velocities  through  the  sluice  gates  and  reduced 
velocities  within  the  sail-through  structure.  An  example  of  the  hydrauli¬ 
cally  efficient  structure  and  the  inefficient  structure  are  shown  in  Figures 
35  and  36,  respectively. 
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Figure  35.  Efficient  sail-through  structure  with  efficient  sluice  gates. 


Figure  36.  Less-efficient  sail-through  structure  with  efficient  sluice  gates. 
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The  initial  testing  phase  was  primarily  intended  to  obtain  the  size  of  the 
structure  and  the  necessary  number  and  sizes  of  the  sluice  gates.  Because 
the  exact  design  was  not  known  at  the  time  of  this  study,  a  general  design 
was  used  in  the  mesh.  A  list  of  the  alternatives  tested  in  the  initial  testing 
phase  is  given  in  Table  l. 


Table  1.  Plan  Alternatives. 


Structure 

Sluice  Gates 

Alternative  Number 

Depth 

(ft) 

Width 

(ft) 

Depth  of 
Sluice  Gate 

1  (ft) 

Width  of 
Sluice  Gate 

1  (ft) 

Depth  of 
Sluice  Gate 

2  (ft) 

Width  of 
Sluice  Gate 

2  (ft) 

Depth  of 
Sluice  Gate 

3  (ft) 

Width  of 
Sluice  Gate 

3  (ft) 

Existing 

Conditions 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

Existing 

Conditions  with 
new  levee 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

1 

12 

56 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

2 

12 

56 

12 

46 

N/A 

N/A 

N/A 

N/A 

3 

10 

56 

12 

46 

12 

46 

N/A 

N/A 

4 

12 

56 

12 

46 

12 

46 

N/A 

N/A 

5 

10 

56 

10 

46 

10 

46 

10 

46 

6 

10 

56 

12 

46 

12 

46 

12 

46 

7 

12 

56 

12 

46 

12 

46 

12 

46 

8 

12 

56 

12 

60 

12 

60 

12 

30 

9 

12 

71 

12 

46 

12 

46 

12 

46 

10 

12 

200 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

2.7  Computational  environment 

The  hydrodynamic  modeling  was  executed  on  the  ERDC  High  Perform¬ 
ance  Computing  (HPC)  SGI  Origin  3000  (Ruby)  parallel  processing  super¬ 
computer.  The  initial  model  mesh  contains  57,834  nodes  and  84,790  ele¬ 
ments,  with  the  mesh  adapting  the  number  of  nodes  to  as  many  as 
140,000  when  needed.  The  model  was  executed  on  16  parallel  processors 
and  required  approximately  20  hours  of  computational  time  (54.6  CPU 
hours)  to  run  for  14  model  days.  The  model  used  a  time  step  of  360  sec¬ 
onds,  with  the  ability  to  adapt  to  smaller  time  steps  when  needed. 
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3  Model  Study  Results 

3.1  Plan  configuration  results 

The  maximum  flood  and  ebb  velocities  in  the  sail-through  structure  are 
plotted  as  a  time  series  in  Figure  37.  The  flow  through  Bush  Canal  is 
driven  by  the  water  level  differences  between  Bayou  Terrebonne  and 
Bayou  Petit  Caillou.  The  smallest  velocities  are  for  the  existing  conditions, 
which  offer  the  best  connections  between  the  two  bayous  and  allow  the 
water  levels  to  equilibrate  most  quickly.  When  the  proposed  levee  is  built, 
many  additional  connections  will  be  cut  off,  making  the  remaining  ones 
pass  more  water  at  higher  velocities.  This  effect  is  exemplified  in  a  model 
run  using  the  existing  conditions  with  only  the  new  levee  represented  in 
the  mesh  (i.e.,  no  structure  in  Bush  Canal).  The  removal  of  the  small  side 
channels  connecting  Bayou  Petit  Caillou  and  Bayou  Terrebonne  (see 
Figure  17)  resulted  in  an  increase  of  0.7  ft/s  in  the  maximum  flood  velocity 
and  a  0.4  ft/s  increase  in  the  maximum  ebb  velocity.  When  this  result  is 
paired  with  a  major  constriction  in  the  channel  (proposed  alternatives),  a 
significant  increase  in  velocity  occurs.  Figure  38  is  a  bar  graph  of  the 
maximum  model  velocities  in  the  structure  for  the  different  alternatives. 


Max  Velocities 
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Figure  37.  Maximum  velocities  in  the  sail-through  structure;  positive  is  flood  direction. 
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Figure  38.  Maximum  model  velocities  in  the  sail-through  structure. 


Model  Velocities  in  the  Sail  Through  Structure 


From  the  initial  model  results,  it  was  determined  that  a  maximum  of  three 
sluice  gates  should  be  used.  As  the  number  of  sluice  gates  increased,  the 
effectiveness  of  the  additional  gates  was  reduced.  The  additional  sluice 
gates  were  farther  from  the  natural  channel,  thereby  reducing  their  effec¬ 
tiveness.  Shown  in  Figure  39  is  a  contour  velocity  plot  of  the  maximum 
flood  conditions  for  alternative  5  (56  ft  sail-through  structure  with  three 
46  ft  sluice  gates).  The  velocities  in  the  third  sluice  gate  are  significantly 
lower  than  the  velocities  in  the  first  and  second  sluice  gates.  From  the  bar 
plot  in  Figure  38,  the  increase  from  adding  the  third  sluice  gate  is  lower 
than  the  increase  gained  from  adding  the  second  gate.  Based  on  this  trend, 
it  may  reasonably  be  concluded  that  the  addition  of  a  fourth  sluice  gate 
would  provide  only  minimal  improvement  in  the  velocity  magnitudes. 
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Figure  39.  Velocity  vectors  and  contours  for  alternative  5. 


3.2  Wind  sensitivity  analysis 

Because  this  model  describes  a  shallow  system,  wind  effects  may  be  im¬ 
portant.  A  sensitivity  analysis  was  performed  to  determine  an  approximate 
effect  on  the  velocities  due  to  wind.  Wind  measurements  for  the  Bayou 
Petit  Caillou  U.S.  Geological  Service  (USGS)  gage,  shown  in  Figure  40, 
were  applied  to  the  numerical  model  to  determine  an  approximate  magni¬ 
tude  of  wind  effects  on  the  currents  through  the  structure.  Because  alter¬ 
native  5  was  judged  to  be  the  best  alternative,  it  was  used  for  this  analysis. 
A  percentile  analysis  was  performed  on  the  model-outputted  velocity  re¬ 
sults,  without  and  with  wind,  for  a  3  month  test  period  from  May  to  July 
2006.  Wind  data  are  plotted  in  Figure  41.  It  was  observed  that  the  maxi¬ 
mum  values  increased  with  wind  by  approximately  15%,  from  -3  ft/s  to 
-3.5  ft/s  (Figure  42  versus  Figure  43).  While  that  increase  is  significant, 
the  90th  percentile  velocities  (i.e.,  velocities  at  or  below  these  values  90% 
of  the  time)  produced  very  few  differences  (Figure  44  versus  Figure  45). 
There  appears  to  be  a  steady  northwestern  wind  blowing  during  this  time 
period  (see  Figure  46).  Because  wind  fields  are  rarely  this  constant,  it  is 
believed  that  the  data  captured  a  rare  event  that  is  unlikely  to  occur  very 
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often.  Therefore,  wind  effects  can  have  a  significant  impact  on  the  veloci¬ 
ties  within  the  sail-through  structure  during  rare  events,  but  they  do  not 
appear  to  produce  a  consistent,  significant  increase  in  the  velocity  fields. 


Figure  40.  Location  of  USGS  gage  for  wind  data. 


Time,  (hrs  from  Jan  1  0:00,  2006) 


Figure  41.  Wind  data  for  the  Bayou  Petit  Caillou  USGS  gage  (USGS  07381343)  for  2006. 
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Figure  43.  Maximum  velocities  for  Alternative  5,  with  wind. 
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Figure  44.  90th  percentile  velocities  for  Alternative  5,  without  wind. 


Figure  45.  90th  percentile  velocities  for  Alternative  5,  with  wind. 
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Figure  46.  Wind  data  during  maximum  velocity  model  results. 
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4  Conclusions 

Using  the  information  developed  in  this  study,  New  Orleans  District  de¬ 
termined  that  alternative  5  produced  the  best  benefit/cost  ratio  with  fa¬ 
vorable  velocity  results.  Alternative  5  consists  of  a  sail-through  structure  of 
56  ft  in  width  and  10  ft  in  depth,  with  three  sluice  gates  each  46  ft  in  width 
and  10  ft  in  depth.  This  alternative  is  illustrated  in  Figure  47.  Other  alter¬ 
natives  produced  better  velocity  results,  but  the  improvement  was  deemed 
not  worth  the  significant  additional  costs  associated  with  those  alterna¬ 
tives.  It  should  be  noted  that  these  velocity  results  are  highly  dependent  on 
the  actual  specifications  of  the  structure  design.  As  discussed  and  illus¬ 
trated  in  Chapter  2,  if  the  sail-through  structure  is  more  efficient  than  the 
one  modeled,  the  velocities  may  be  higher  than  those  reported  here  for  the 
same  structure  combination.  Another  important  consideration  is  the  effect 
of  wind,  which  can  be  significant,  but  is  not  believed  to  impact  the  opera¬ 
tion  of  the  structure  represented  by  alternative  5  in  normal  conditions. 
Significant  wind  events  large  enough  to  notably  change  the  velocities  are 
believed  to  be  a  rare  occurrence. 


Figure  47.  Alternative  5. 
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Appendix:  Description  of  the  ADaptive 
Hydraulics  (ADH)  Model 

ADH  is  a  state-of-the-art  ADaptive  Hydraulics  Modeling  system  developed 
by  the  U.S.  Army  Engineer  Research  and  Development  Center  -  Coastal 
and  Hydraulics  Laboratory  (ERDC-CHL).  It  is  capable  of  simulating  both 
saturated  and  unsaturated  groundwater,  overland  flow,  three-dimensional 
Navier-Stokes  flow,  and  two-  or  three-dimensional  shallow  water  prob¬ 
lems  with  the  current  study  utilizing  the  two-dimensional  shallow  water 
module.  The  2D  shallow-water  equations  used  for  this  application  are  a 
result  of  the  vertical  integration  of  the  equations  of  mass  and  momentum 
conservation  for  incompressible  flow  under  the  hydrostatic  pressure  as¬ 
sumption.  Written  in  conservative  form,  the  2D  shallow  water  equations 
are: 


5U  5F  SG 

—  +  —  + - + 

dt  8x  dy 


H  =  0 


where 


\h  1 


u  = 


uh 


w 


(1) 


(2) 


uh 

F  =  \u2h  +  \gh 2  - 
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uvh  -  h  — — 
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vh 
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and 


where: 


H 


0 


,  dz,  2  ,  «V«2  +  v2 
gh  J  +  n2gh  3 

ox  C  0h 

,  dzt  ,  2  ,  vVi 

gh  — -  +  «  gh - - 

sy  cy 


Vy/ll2  +  V2 


(5) 


p  =  fluid  density 

g  =gravitational  acceleration 

Zb  =  riverbed  elevation 

n  =  Manning's  roughness  coefficient 

h  =  flow  depth 

u  =  x-component  of  velocity 

v  =  y-component  of  velocity 

Co  =  dimensional  conversion  coefficient  (l  for  SI  units,  1.486  for  U.S. 
customary  units) 

a's  =  the  Reynolds  stresses  due  to  turbulence,  where  the  first  subscript 
indicates  the  direction,  and  the  second  indicates  the  face  on 
which  the  stress  acts. 


The  Reynolds  stresses  are  determined  using  the  Boussinesq  approach  to 
the  gradient  in  the  mean  currents: 


a 


XX 


2  pv 


du 
1  dx 


avy  =  2PV 


dv 
'  dy 


and 


^.vv  =  =  2  pv, 


f  du  dv  ^ 
dy  dx 


(6) 

(7) 


(8) 


where  vt  =  kinematic  eddy  viscosity  (which  varies  spatially). 
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The  ADH  shallow- water  equations  are  placed  in  conservative  form  so  that 
mass  balance  and  the  balance  of  momentum  and  pressure  are  identical 
across  an  interface.  This  is  important  in  order  to  match  the  speed  and 
height  of  a  surge  or  hydraulic  jump. 

The  equations  are  represented  in  a  finite  element  approach.  The  quality  of 
the  numerical  solution  depends  on  the  choice  of  the  basis/trial  function 
and  the  test  function.  The  trial  function  determines  how  the  variables  are 
represented  and  the  test  function  determines  the  manner  in  which  the  dif¬ 
ferential  equation  is  enforced.  In  the  Galerkin  approach  the  test  functions 
are  chosen  to  be  identical  with  the  trial  functions.  When  the  flow  is  advec- 
tion-dominated,  the  Galerkin  approach  produces  oscillatory  behavior.  The 
Galerkin  form  of  the  test  function  cannot  detect  the  presence  of  a  node-to- 
node  oscillation  and  so  allows  this  spurious  solution.  The  approach  used  in 
ADH  is  to  enrich  the  standard  Galerkin  test  function  with  an  additional 
term  that  can  detect  and  control  this  spurious  solution. 


This  Petrov-Galerkin  method  used  here  is  based  on  elemental  constants 
for  coefficients.  This  reduces  the  stabilization  to  the  nonconservative  form. 
This  is  not  a  problem  for  conservation  since  the  stabilization  is  only  ap¬ 
plied  within  the  elements  and  uses  the  Galerkin  test  function  to  enforce 
“flux”  balance  across  element  edges.  For  illustration,  consider  the  shallow- 
water  equations  in  nonconservative  form 


dU  AdU 

- +  A - 

dt  dx 


dy 


+  H 


0 


(9) 


where  A  =  —  and  B  =  — .  The  trial  functions  (or  interpolation/basis 

dU  dU 

functions)  are  the  Lagrange  polynomials.  These  are  piecewise  linear  func¬ 
tions  that  are  continuous  across  element  boundaries.  Spatial  derivatives, 
however,  are  not  continuous  across  these  element  edges.  Each  of  the  de¬ 
pendent  and  independent  variables  is  interpolated  via  these  trial  func¬ 
tions.  For  example, 


N 

u(x)  =  YJ<Pj{x)uj 
j= i 


(10) 
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means  that  the  approximate  solution  is  made  up  of  the  sum  of  the  prod¬ 
ucts  of  the  trial  function  for  node  j  and  the  nodal  value  at  that  location. 
The  test  function  is  chosen  as: 


where, 


</>*=</>,!  + a 


dx 


dy 


\ 

J 


r  /  > 

2 

a  =  0.5/ 

v  •  v  +  gh  + 

V  At  J 

and 

1  =  (fie)1/2,  the  square  of  the  element  area 
v  =  ( u  ,  v ),  the  element  average  velocity  components 
At  =  time  step  size. 

The  finite  element  statement  becomes: 


(ll) 


(12) 


(13) 


where  the  subscript  /  indicates  the  finite  element  approximation.  The  Pet- 
rov-Galerkin  contributions  are  integrated  over  the  interior  of  the  elements, 
but  not  across  element  edges.  This  contribution  stabilizes  the  Galerkin  ap¬ 
proach.  This  scheme  utilizes  a  single  scaling  factor,  a  .  This  is  different 
from  the  scheme  reported  in  Berger  and  Stockstill  (1995).  That  scheme  in¬ 
volved  scaling  each  eigenvalue,  but  that  method  does  not  converge  using 
the  iterative  solver  in  ADH.  Instead,  a  single  value  scaling  (Equation  12)  is 
used. 


One  of  the  major  benefits  of  ADH  is  its  ability  to  adapt  the  mesh  in  areas 
where  additional  resolution  is  needed  to  properly  resolve  the  hydrodynam- 
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ics.  This  process  is  done  by  normalizing  the  results  so  that  an  error  quan¬ 
tity  is  determined  for  each  element.  If  this  error  exceeds  the  tolerance  set 
by  the  user,  then  the  element  is  refined.  ADH  is  also  able  to  unrefined  pre¬ 
viously  refined  areas  when  the  added  resolution  is  no  longer  needed.  ADH 
contains  other  essential  features  such  as  wetting  and  drying,  completely 
coupled  sediment  and  salt  transport,  and  wind  effects.  A  series  of  modu¬ 
larized  libraries  make  it  possible  for  ADH  to  include  vessel  movement, 
friction  descriptions,  as  well  as  a  host  of  other  crucial  features.  ADH  can 
run  in  parallel  or  on  a  single  processor  and  runs  on  both  Windows  systems 
and  UNIX  based  systems. 
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